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Introduction 



Recently, much attention has been paid to the study of integrable field theory with a 
boundary The study of such field theories is not only intrinsically 

interesting but also provides a better understanding of boundary related phenomena in 
statistical physics and condensed matter || [|TIJ . Probably the most famous physical ex- 
ample of a boundary integrable model is the Kondo problem, where the 1+1 dimensional 
field theory is an effective field theory of s-wave scattering of electrons off a magnetic 
spin impurity. Such an impurity problem, in which one concentrates on s-wave scattering 
from some isolated object at the origin, generically provides interesting 1+1 dimensional 
boundary field theory. There are also 1+1 dimensional boundary quantum integrable 
systems of experimental relevance, such as a Luttinger liquid (Thirring) model for the 
edge states of electrons in the fractional quantum Hall effect |i0| . 



An integrable field theory possesses an infinite set of independent, commuting inte- 
grals of motion. In the 'bulk theory' these integrals of motion follow from an infinite 
number of divergenceless currents. However, when the theory is restricted to a half- 
line (or to an interval) the existence of conserved charges on the whole-line does not 
guarantee integrability unless special boundary conditions are specified. The boundary 
conditions under which a theory preserves its integrability can be sought in several ways: 
for example, via a perturbed conformal boundary condition [ njlHl2 | or via a Lax pair 
approach [|J. 



However, it is more difficult to discuss the integrability for non-ultralocal integrable 
models, such as a nonlinear sigma model on the half-line, because their fundamental 
Poisson bracket relations are field dependent |L3| ]14[ |15| Jl6|] [[L7| |18 . Ghoshal[19j has 
obtained the boundary S-matrix of the 0(N)-symmetric nonlinear sigma model, when 
he conjectured that this model is still integrable in the semi-infinite space with the 
'free' or 'fixed' boundary conditions. However, Mourad and Sasaki[|2l| have concluded 
recently that the infinite set of non-local charges characterizing integrability on the whole 
line is not preserved even for the free (Neumann) boundary condition. In this paper, 
the question of integrability for sigma models with a boundary will be re-examined 
from a different view-point. First, standard tools for the model on the full line will 
be reviewed. In particular, the Lax pair representation, the r, s matrices, and their 
use to generate conserved charges, will be discussed in the next section. Following the 
ideas of [f| , a Lax pair on the half-line will be constructed in the third section, with the 
property that its zero curvature condition leads to both the equation of motion and to 
the boundary condition. In the fourth section, the Lax pair will be used to construct 
the conserved charges, suggesting integrability of the model on the half-line, and to 
determine compatible integrable boundary conditions. 



2 



The 0(N) nonlinear sigma model on the full line 



To establish notation, the standard Lax pair for the SO(N) nonlinear sigma model 
13| will be written in the form 



a = T A?(\A + A 1 ) 
ax = T ^{XA 1 + Aq) 



where 



A ll . = 2(n®d ll n-d ll .n®n), (2) 



and n is an iV-component, unit length vector. 



It is easy to check that the zero curvature condition for the vector potentials defined 
by eq(l), 

d ai - dia + [ax, a ] = 0, (3) 
leads to the equation of motion of the model in the form 

n (g> d 2 n — d 2 n ® n = 0, (4) 

or, equivalently, 

d 2 n a + ri a ((9 M n 6 <9 M n 6 ) = 0. (5) 
The action of the O(N) nonlinear sigma model may be taken to be 

S = \jd 2 x {{d,n a f + CK - 1)) , (6) 

where ( is a Lagrange multiplier. The canonical momentum associated to n a is 7c a = don a . 
However, the Hamiltonian formalism is of constrained type. The Hamiltonian is 

H=\J dxK + {d ina f} (7) 

and the canonical Poisson brackets are 

{n a (x), n b {y)} = 

{n a (x), n b (y)} = (5 ab - n a n b )5(x - y) (8) 
{7r a (x), n b (y)} = (n a n b - n b n a )S(x - y) 

It is also useful to introduce a representation for the SO(N) generators 

T$ = S a X - 8%5l (9) 
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with quadratic Casimir 
and to set, using the components (2), 



C = -T ab g> T a \ (10) 



A fM = 2n a d fl n b T ab . (11) 
It is also convenient to define an additional quantity 

A = n a T ab <g> T bc n c . (12) 

The two operators C and A play a significant role in setting up the Poisson structure of 
the model. The use of the same notation for the matrices, or components should cause 
no confusion. 



Using the Poisson brackets (8), and following closely Bordemann et al. [[UJ, one may 
obtain the following algebraic structure: 

{A (x) f A (y)} = 2[C, I ® A )5(x - y) 

{A 1 (x)°A 1 (y)} = (13) 
{A (x) f A x (y)} = 2[C, I ® A^ix - y) - AA(y)5'(x - y). 

For the last of these, use has been made of the identity 

(f(x)-f(y))6'(x-y) = -f'(x)5(x-y). 

In turn, eqs(13) may be used to calculate the Poisson brackets of the Lax pair oper- 
ators themselves, casting them into the Maillet form: 

{ai(x, A) f ai(y, fi)} = {[r(x, A, fi) , ai(x, A) <g) I + I ® ai(x, /x)] 

— [s(x, A, /i) , ai(x, A) (8) J — I ® ai(x, /i)])<5(x — y) 
+ (r(x, A, /i) + s(x, A, /i) - r(y, A, /i) + s(y, A, - (14) 

where 

r{x X a)- 2A ^ C + 2A/i(A-/i)(l + A//) 

1 ' ' *° ~ (A - - A/i) 6 + (1 - - /i 2 )(l - A/*) ( } (15) 

s(x, A, /i) = j^z^n^-^ A i x )- ( 16 ) 
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In checking eqs(14) the following additional relationship is indispensable 

[A, I®A li -A li ®I] = [C, Ap®!] (17) 



In equations (14), r(x, A, /i) and s(x, A,/i) are two matrix 'structure constants' which 
depend on x because they depend on the field of the theory. They have following property 



r(x, A, fi) = —r(x, —A — /i) 
and satisfy the following modified Yang-Baxter equation 

(r + s)2a(x,fi,ri) , (r + s) 12 (x, A, 77) + (r + s) 23 (z, H, v) , (r + s) 13 (x, A, 77) 



(18) 



+ 



(r + s) 13 (a;,A,77) , (r - s) 12 (x, A,//) + iff^x, A, n, 77) - H^ 3 s> (x, A, 77) = 0, 



r(H-s), 



(19) 



where ifj^f (2;, A, //, 77) is defined by 

{ai(x,A)(g)7 , (r + s) 23 (y, /X, 77)} = H[^{x,X,iJL,r])S(x - y). (20) 



One of the attractive features of the Lax pair representation is the role it plays in 
the generation of conserved quantities. The path-ordered exponential 



(21) 



(22) 



rv 

T(x, y; A) = P exp / ai(z,X)dz 

Jx 

enjoys the following properties, 

T(x,x; A) = 1, T(y, x; A) = T~ 1 (x, y; A), 
T(x, y; X)T(y, z; A) = T(x, z; A), 

and satisfies 
{T(x, y\ A) f T(x', y'; //)} 

= /» cfe J*' <fc'T(x, z- A) <g> T(x', z'; fj,){ ai (z, A) ® a^', /i)}T(z, y, A) <g> T(z', y'; //). 

(23) 

Defining £/(A) = T(— 00, +00; A), one obtains from eq(23) the relation 

{[/(A) ® [/(//)} = r(-cx), A, n)U(X) ® U(/i) - U(X) ® U(fi)r (+00, A, //). (24) 

Alternatively, assuming periodic boundary conditions on the finite interval (—L,L), ie 
n a (—L) = n a (+L), one has 



{U L (\)®U L (n)} = [r(-L,\,») , C/l(A)®C/ l (//)] 
+(C/l(A) ® I)s(L, A, /i)(J ® UM) -(I® U L (fi))s(L, A, /i)(t4(A) ® I), 



(25) 



where Ul(X) = T(—L, L, A). From eq(24) or eq(25), it is straightforward to obtain 

{trU(X) ,trU{n)} = (26) 

indicating that the infinity of conserved charges obtained from U (A) as coefficients in an 
expansion in powers of A are in involution. 



Lax pair representation for O(N) sigma model on a half-line 

A Lax pair on a half-line will be constructed so that its zero curvature condition 
leads both to the field equation and to the boundary condition at x = a. The procedure 
developed in Q will be followed closely. To construct a modified Lax pair including the 
boundary condition, it is first of all convenient to consider an additional special point 
x = b (> a) and two overlapping regions R_ : x < (a+5+e)/2; and R + : x > (a+b—e)/2. 
The second region will be regarded as a reflection of the first, in the sense that if x G R + , 
then 

n a (x) = n a (a + b - x). (27) 

The regions overlap in a small interval surrounding the midpoint of [a, b]. In the two 
regions, define: 

: do 
a i 

di 

It is then straightforward to check this Lax pair leads to the equation of motion on the 
half-plane, and to the boundary conditions at the boundary point x = a: 

d x n a = -T a . (29) 

For the moment, the question of the form of the boundary potential term which ought 
to be added to the Lagrangian will be postponed. 

On the other hand, for x G R- and x > a, a\ vanishes and therefore the zero curvature 
condition merely implies that do is independent of x, which means in turn that the fields 
n a are independent of x in this region. Similar remarks apply to the region x G R + 
and x < b. Hence, taking into account the reflection principle, fields are independent 
of x throughout the interval [a, b], and equal to the field value at a or b. For general 
boundary conditions, the gauge potential do is different in the two regions R±. However, 



ao + x^i) 6 ^ x (din + J 7 ) - (d x n + J 7 )<g>n] 

9 (a — x)a\ 

ao + 1^)0(0 - x) [n ® (d x n - T) - (d x n -F)®n\ 
9(x — b)a x 



(28) 
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to maintain the zero curvature condition over the whole line the values of So must be 
related by a gauge transformation within the overlap region [a, b]. In principle, this gauge 
transformation could be quite general, even depending upon the field n. Here, it will be 
assumed a gauge transformation k(A) exists, not depending upon either n, or x°. Thus, 

8qk = nao(x°, b) — do(x°, cl)k = 0. (30) 

These are strong assumptions. Previous experience with affine Toda theory Q indicates 
they might not lead to stronger constraints on the boundary condition than those which 
would be obtained by alternative means (such as examining directly the existence of 
conserved quantities). However, the classical r-matrix for affine Toda field theory does 
not itself depend on the Toda field. This fact represents a crucial difference. Substituting 
do into the above, one obtains 

X[k(X),A ]_ + [k(X),B} + = (31) 

where 

B = 2(n®F-F®n). 

Obviously, k(X) = I and T = is one of the solutions; it corresponds to the free boundary 
condition. To find the others for general N is less easy and only the case N = 3 will be 
discussed in detail here. 



It is useful first to note that when N = 3 k(X) is an element of the group 50(3) and 
must therefore have the form 

Kab = PSab + QKh + Re a b C k c = Kab + K ib ( 32 ) 

where P, Q, R and the vector k could depend upon A but not n or x°, and the S, A 
superscripts denote the symmetric and antisymmetric parts, respectively. Orthogonality 
requires 

p 2 + \k\ 2 R 2 = 1 \k\ 2 Q 2 - R 2 + 2PQ = 0. (33) 
It is also useful to write 

(Ao)ab — tabc^c B = e a b c (3 c . 

Taking the symmetric and antisymmetric parts of eq(31) yields a pair of equations: 

* ( K ad e db C a c + K-M e da C a C ) + R (Pah + PbK) = 

XR (a a k b - a b k a ) + n s ad e dbc f3 c - ^ d e dac (3 c = 0. 

The vanishing trace of the first of eqs(34) requires R (/3 • k) — 0. Therefore, assuming 
R 7^ 0, (3 must have the form 

(3 = pxk. (35) 
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The second of eqs(34) forces p to be proportional to a and, without any loss of generality, 
(3 may be set equal to a x k. However, it is also necessary that 

XR + 2P+ \k\ 2 Q = 0. (36) 

Finally, returning to the first of eqs(34) and using (35) the other components vanish 
provided 

\Q + R = 0. (37) 
Eqs(33), (36) and (37) (over-) determine P,Q and R to be 

P - ± A^T^ AH^F AM^' (38) 

Finally, noting 

(3 = 2nx T a = 2nx d n, (39) 

and, therefore, 

n x T = (n x <9 n) x fc, (40) 

it is straightforward to deduce 

T = -k x <9 n + (n ■ k x 9 n)n (41) 
together with the further constraint: 

k-d n = 0. (42) 

To summarise, a consistent boundary condition requires 

din = — k x <9 n + (n • k x d n)n and • <9 n = at x = a. (43) 

The vector may be chosen freely and represents extra parameters allowed by the 
boundary condition. Notice, however, that with the exception of the Neumann condition, 
the boundary violates the 0(3) symmetry of the model leaving a subgroup of rotations 
preserving k. 

The choice R = leads to T = and k — I. 

Construction of conserved quantities 

The patching matrix k(X) is an essential ingredient in the construction of conserved 
quantities for the model defined on the half-line x < a via the trace of the generating 
function: 

U(X) = T(-oo, a; A)/c(A)T(6, oo; A). (44) 
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Using the reflection properties, 

n(x) = n(a + b — x) d\n(x) = —din(a + b — x) 

d n(x) = d n(a + b — x) ai(x, A) = —a^a + b — x, —A), 
eq(44) may be rewritten as 

U(X) = T(-oo, a; X)k(X)T-\-oo, a; -A). (45) 

Integrability requires the conserved quantities to be in involution and, when the 
model is restricted to a half-line, this in turn implies a compatibility relation which must 
be satisfied. This relation is expected to involve k(X) and the classical r matrix entering 
eq(14). The Poisson brackets: 

{T(x, y; A) f T-\x, y; //)} = (T(x, y; A) <g> /)r(y, A, //) (/ ® T _1 (x, y; »)) 

- (/ <g) T- 1 (x, y; /i)) r(x, A, /x) y; A) ® i) , (46) 

{T- 1 ^, y; A) ® T(x, y; //)} = -{T~\x, y; A) ® l)r(x, A, /i) (/ ® T(x, y; A)) 

+ (/ ® T(x, y; A))r(y, A, /i) (t- 1 (rr, y; //) ® /) , (47) 

and 

{T^^yjA) ® T- 1 (x,y;/i)}= (T- 1 (x,y;A)®T~ 1 (x,y;/i))r(a;,A, / u) 

- r(y, A, n) (T-\x, y; A) ® T" 1 (x, y; /i) ) (48) 

imply the following Poisson bracket for a pair of generating functions, as defined in eq(45) 
but with different parameters A and /i, 

MA) f [/(//)} = r(-oo, A, /i) (C/(A) ® [/(//)) - (u(X) ® C/(//))r(-oo, A, //) 
+ ([/(A) ® /)r(-oo, A, -/i)(/ ® [/(a*)) - (/ ® f/( / u))r(-oo, A, -/i)(f/(A) ® /) 
+ (t(X) ® T(fi)^ —r(a, A, fi) k(X) ® /t(/i) + «(A) <g> «;(//) r(a, A, //) 

+ (/ ® /t(/i))r(a, A, -/i) (k(A) ® /) - (k(A) ® /)r(a, A, -//) (/ ® re(^)) • 

.(T-^-^OT-H-Zi)) 

(49) 

where T(A) = T(-oo, a; A). To establish 

{trU(X) , trf/(/i)} = (50) 
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(51) 



it would be sufficient to require 

k(X) (g> «;(//), r(a, A,//) + (i ® «(//)Jr(o, A, — ^)(«(A) <8> /J 

- (/c(A) ® l)r(a, A, -/*)(/ = 0, 

a condition reminiscent of those encountered in [JTJ] and ||. Obviously, eq(51) is satisfied 
by k(A) = /, corresponding to the free boundary condition. However, it is not in general 
satisfied for the other possibilities, principally because k is independent of n while r 
depends explicitly on the field n but not its derivatives. 



Alternatively, to demonstrate eq(50) it would be sufficient to require 



k(X), T-\-\)T(\) 



0. 



However, while no proof (or disproof) of this relation is available, it too would appear 
to be unlikely to be true given the conjectured n-independence of k. 



Conclusion and Discussion 



It has been argued that at the classical level the O(N) nonlinear sigma models remain 
integrable when restricted to a half-line with a Neumann boundary condition. Infinitely 
many conserved charges in involution have been constructed via eq(41) and one might 
suppose this result to be a strong indicator of integrability. However, Mourad and 



Sasaki |2T[ have taken a different view, analysing the non-local charges for the theory 
restricted to a half-line, and concluding the opposite. Their discovery of the absence 
of non-local charges in the presence of the Neumann condition would appear to supply 
strong reasons for a lack of integrability for even the simplest of boundary conditions. 
On the other hand, others (for example, Ghoshal QT9J and Mackay f20| ) have assumed 



the sigma models restricted to a half-line are quantum integrable, and have conjectured 
quantum reflection matrices compatible with the reflection bootstrap. One would expect 
the quantum integrability to be reflected in classical features of these theories (although 
it must be borne in mind that the quantum sigma models have a spectrum and other 
properties which are very far from being reflected in the classical degrees of freedom). It 
would seem one might be forced to conclude that the classical non-local charges are not 
crucial for integrability. 



It has been shown that the 0(3) sigma model on a half-line has a consistent Lax pair 
description not only for the Neumann condition at the boundary, but also for some other 
suitably chosen conditions, reported in (43). However, it is not yet clear if these alter- 
native boundary conditions allow charges in involution. These conditions also present 
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something of a puzzle because they do not appear to be derivable from a boundary po- 
tential added to the bulk Lagrangian. The condition (41) would be fine by itself and 
follow from the boundary potential 

B= -n- k x d n+ Un 2 - 1), (52) 

where r is a Lagrange multiplier. However, the additional constraint (42) needs to be 
imposed 'by hand'. The significance of this is not really clear. However, this fact and 
the remarks above suggest that the Neuamnn condition is the only boundary condition 
compatible with the classical Poisson structure. 
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